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FOREWORD 

This  technical  report  documents  a  computer  algorithm  for  the  "influence  num¬ 
bers'  that  govern  the  relationship  between  certain  subsets  ^  rbe  earth's  gravita¬ 
tional  potential  expansion  coefficients  and  "lumped  coefficients,"  each  of  which 
summarizes  the  gravitational  effect  of  an  entire  group  of  potential  expansion 
coefficients  on  the  inclination  resonance  of  earth  satellites. 

The  work  was  performed  in  the  Space  and  Surface  Systems  Division  as  part  of  a 
larger  effort,  with  the  purpose  of  establishing  equations  of  condition  for  groups 
of  gravitational  field  parameters  as  a  contribution  to  the  validation  of  the 
WGS-84  Geodetic  Parameter  Solution. 


The  present  report  was  reviewed  by  C.  W.  Duke,  Jr.,  Head  of  the  Space  and 
Surface  Systems  Division;  R.  L.  Kulp,  Head  of  the  Space  and  Ocean  Geodesy  Branch; 
and  Dr.  C.  Oesterwint er ,  Research  Associate. 
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INTRODUCTION 

An  artificial  earth  satellite  is  said  to  be  in  \i/  a  resonance*  with  the 
earth's  gravitational  field  when  it  performs  I?  nodal  periods  while  the  earth  re¬ 
volves  a  times  relative  to  the  satellite's  precessing  orbit  plane.  Near  reso¬ 
nance,  the  perturbations  caused  by  those  tesseral  components  of  the  earth's  gravi¬ 
tational  potential  with  which  the  satellite  happens  to  be  in  resonance  generally 
exceed  in  magnitude  all  the  remaining  tesseral  perturbations. 

Near  resonance,  the  rate  of  change  of  the  orbital  inclination  of  the  satel¬ 
lite  may  be  expressed  as  an  aggregate  of  several  more  or  less  complicated  func¬ 
tions  involving  the  parameters  characterizing  the  resonance  type,  the  orbital 
elements,  certain  quantities  named  "lumped  coefficients,"  and  time.  The  lumped 
coefficients  are  linear  expressions  consisting  of  the  scalar  products  of  sets  of 
expansion  coefficients  of  the  earth's  gravitational  potential  (in  the  sense  that 
each  lumped  coefficient  is  associated  with  a  specific  potential  term  order)  times 
a  second  set  of  parameters  customarily  called  "influence  numbers."  The  latter  are 
essentially  functions  of  orbital  semimajor  axis,  inclination,  and  eccentricity. 

Satellites  with  sufficiently  short  orbital  periods  will,  as  their  orbits 
decay,  pass  through  resonance  or,  frequently,  through  a  succession  of  different 
resonances.  While  orbital  resonance  effects  are  not  restricted  to  the  orbit 
inclination,  conditions  for  acquiring  useful  inclination  resonance  data  have  been 
especially  favorable,  giving  rise  to  a  large  body  of  investigations  (References  1 
through  78).  Most  of  these  are  undertaking  to  evaluate  the  lumped  coefficients 
for  empirical  time  sequences  of  resonant  inclination  values.  There  is  almost 
always  an  attempt  made  to  subsequently  infer  numerical  values  for  the  geopotential 
expansion  coefficients  related  to  the  lumped  coefficients  found.  Most  of  the 
investigations  seen  are  confined  to  specific  resonance  cases  and  to  specific  satel¬ 
lites.  Some  however  attempt  a  synopsis  (see  especially  Reference  (■>/).  An  effort 
was  recently  made  at  t  he  Naval  Surface  Weapons  Center  (NSWC)  to  <-v.st  >«mat  i  r.a  I  1  y 
review  the  entire  body  of  published  research  work  and  to  compile  all  lumped 
coefficients  listed.  A  system  of  linear  conditional  equations  complete  with  the 


customary  statistical  information  was  generated  lor  the  geopot ent.  i  a 
*For  detail...  -ee  tin  li  n  sv  oi  resonance  in  Re!  e  fences  md  .«  ’> . 
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coefficients,  the  latter  being  the  system  unknowns,  and  ranging  up  to  degree  42, 
each  identified  lumped  coefficient  contributing  one  equation,  A  general,  numeri¬ 
cal  solution  for  those  geopotential  parameters  contributing  to  the  inclination 
resonance  cases  treated  in  the  just-mentioned  collection  of  literature  will  in  due 
course  be  executed. 

As  already  stated,  a  necessary  ingredient  of  the  conditional  equations  are 
the  influence  numbers.  These  are,  essentially,  geometry  factors  reflecting  the 
orbital  parameters  at  resonance.  As  a  generally  valid  theory  of  the  influence 
numbers  could  not  be  located  in  any  of  our  sources,  it  was  derived  starting  from 
the  Lagrange  planetary  equation  for  the  inclination  variation.  To  document  the 
resulting  computer  algorithm  is  the  purpose  of  the  present  report. 
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ALGORITHM  FOR  THF  INFLUENCE  NUMBERS 


INPUT  DATA 


R  .  "Radius  of  the  Earth"  (semimajor  axis  of  reference  ellipsoid) 

in  km.  Default  value:  R  =  6378.140  km. 

a  .  Semimajor  axis  of  satellite,  in  km. 

e  .  Orbital  eccentricity 

i  .  Orbital  inclination 

3  1  Integers,  characterizing  3/a  orbital  resonance, 

a  1  .  Normally,  a  =  1  or  2.  In  most  cases,  3  will  be  larger  than 


5  und  smaller  than  50. 

Integer,  denoting  magnitude  of  the  resonance  term  (order  of 
"overtone").  Generally,  y  will  be  equal  to  1  or  2.  Generally, 
it  will  be  a  positive  integer. 

Integer  (positive,  negative  or  zero)  associated  with  the  magnitude 
of  the  orbital  eccentricity.  For  the  present  purpose,  the 
following  values  of  q  may  occur:  0,  ±1 ,  ±2,  ±3. 


q 
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THE  INFLUENCE  NUMBERS 

k  =  ay  -  q 


m  =  3y 


K=  (S  -  c*)y  +  q 

If  K  is  even 


i  =  By  +  2r  =  m  +  2r 


m  =  3y 
r 


+ 


r 


r 


If  K  is  odd 


t'  =  By  +  (2r  +  1 )  =  tn  +  1  +  2r 
r 


Pr  =  ■y(3  -  ct)y  +  |  q  +  7  +  r 


=  -;-(K  +  1)  +  r 


r  -  0, 


r 

max 


i'/.peeti  d  maximum  value  for  r 

1  m;i  y 


r  <  50 
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The  influence  numbers  are  then 


(-Dr(Vr  ^ 

Q  (.*  i, 


FAN 


Imp 


»  FAN, 

Voq  Fompo 


where 


t  =  a 

r 

p  =  p 
i  r 

eQ  =  t(r  =  0) 

PQ  =  P(r  =  0) 

Note  that  the  influence  numbers  are  dimensionless  quantities 


INCLINATION  FUNCTION  SUBROUTINE 


FA 


imp 


=  FAfmp(i) 


a  +  m)  ! 
2fp!(f  -  p)! 


max 

E 

n  . 

min 


*  c3l'-m-2p-2r,  *  sm-('+2p+2r; 
,i  .. 

L  =  cos  (—  l  ) 

.  J 

S  =  S  l  II  ( -  I  ) 

•  .  =  rn.ixj  0,  ()'-m-2p)  | 


-  nun  [  (l-n)  ,  (2  1-2  p)  1 


(13) 

(iA) 

,17) 

i  i  -  >  >!• 
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Stop  summation  over  k  when  [G],,., 


and  [G]XI  agree  to  six  significant  figures 
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(***■*<*) 


For  q'  >  0,  h  =  k  +  q' 


For  q  <  0 ,  h  =  k 


w,,  .  e  y 

/  n  n  ^ 


( "2p'  \  J_  /  U'  -  2p/  +  q')e\r 

\h  -  7  r!  \  2(5  / 


For  q'  >  0,  h  =  k 


For  q  <-  0,  h  =  k  -  q' 


To  avoid  difficulties  arising  from  use  of  negative  arguments,  calculate  the 
binomial  coefficients  occurring  in  Equations  26  and  29  as  follows: 


(b)  *  Bi 


I.  A  -  0 1  *  [  A  -  1  ]  *  [  A  -  2  ]  *  .  *[A  -  (B  -  1)1 


B  may  be  expected  to  be  a  positive  integer  always.  A  may  assume  negative  values. 


ALGORITHM  FOR  I’HE  LUMPED  COEFFICIENTS 


Lot  C  and  S  .  bo  the  normal i/ect  expansion  coo  I  i  i c i on t s  of  the  goonotenti. 
t  m  t  m  1 


and  C  and  S  their  umvrmal  i/e-.l  counterparts,  then 

■.  m  ■  iv. 
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where  N^m  is  the  normalization  factor  specified  bv  Equation  19.  The  relationship 

~*—g  ^  =n  ic 

among  the  lumped  coefficients,  C^’  and  S^’  ,  the  potential  expansion  coefficients 
and  the  influence  numbers  is 
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Note  once  more  that  the  influence  numbers  as  well  as  the  potential  expansion 
coefficients  are  dimensionless. 
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The  following  derivation  should  be  viewed  in  context  with  tin.  i  r.t  roduct  i  <  n 
fo  this  report .  Lagrange’s  planetary  equation  for  the  tnel  in.-tt  i<.»  \  a  r  i .  1 1  i  -mi  ’.;;i 
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..here  a,  e,  i,  „  ,  are  the  familiar  orbital  elements,  M  is  the  mean  motion,  R  is 

flu:  equatorial  "radius  of  Lite  earth,"  M  is  the  mean  anomal;-,  S  is  the  side  real 

Greenwich  time  of  the  epoch  of  the  elements,  and  ('  and  S„  are  the  utitmtmal  >7t-d 
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expansion  coefficients  of  the  gravitational  potential. 

Also  from  Page  288  of  Reference  33,  for  /  a  resonance  tiler.;-  held  tin  relationships 
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as  well  as  Equations  3  through  11  from  the  main  body  of  this  report.  For  the 
meaning  of  beta,  alpha,  gamma,  and  q  see  the  segment  "input  Pat  a"  ire:'.  the 
"Algorithm  for  the  Influence  Numbers. "  Also  note  that 
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If  one  now  enters  Equations  B 1 8  and  B20  respective  Lv  ini  o  Equat.ii  n  Hi  3  ana 
sequent  lv  Equation  Bid  into  l.agi  ungc  '  .->  pianeL.a'v  <.  quatiaii  ( 15<  <  i ,  .....  erv  ny 
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